Electron behavior in topological insulator based P-N overlayer interfaces 
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Topological insulators (TIs) are novel materials that manifest spin-polarized Dirac states on their 
surfaces or at interfaces made with conventional matter. We have measured the electron kinetics of 
bulk doped TI Bi 2 Se3 with angle resolved photoemission spectroscopy while depositing cathodic and 
anodic adatoms on the TI surfaces to add charge carriers of the opposite sign from bulk dopants. 
These P-N overlayer interfaces create Dirac point transport regimes and larger interface potentials 
than previous N-N type surface deposition studies, revealing unconventional Rashba-like and surface- 
bulk electron interactions, and an unusual characteristic distribution of spectral weight near the 
Dirac point in TI Dirac point interfaces. The electronic structures of P-N doped topological interfaces 
observed in these experiments are an important step towards the understanding of solid interfaces 
with topological materials. 



Bismuth based topological insulators (TIs) are the first 
experimentally realized three dimensional topologically 
ordered bulk solids, and feature robust, massless spin- 
helical conducting surface states that appear at any in- 
terface between a TI and normal matter that lacks the 
topological order [H-Q- This topologically defined sur- 
face environment has been theoretically identified as a 
promising platform for observing a wide range of new 
physical phenomena, and possesses ideal properties for 
advanced electronics such as spin-polarized conductivity 
and suppressed scattering jj4l8| . Developments from the 
last year have demonstrated that it is possible to fabri- 
cate nanodevices using_Bi 2 Se3 and pass current through 
the surface states (l9l l20j, but have also highlighted the 
fact that topological insulators are a new state of mat- 
ter distinct from normal metals and insulators, and their 
properties in a junction /interface setting are not known 
or understood 
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In this Letter, we explore this critically important 
question by using cathode surface overlayers with weak 
and strong spin-orbit interactions (Cu and In) and anodic 
nitrogen dioxide in tandem with control of bulk doping in 
TI Bi2Se3 from P-type to N-type charge carrier regimes 
for the first time to generate a range of topological insu- 
lator P-N interface scenarios that are relevant for device 
development 0, H H3- The P-N scenario allows us 
to access the novel aspects of the physics not observed 
in previous studies limited to N-N scenarios. Our ex- 
perimental results reveal diode-like configurations that 
can manifest a gap in the interface electron density near 
a topological Dirac point and systematically modify the 
interface Dirac velocity, allowing far reaching control of 
helical electrons at the interface. In addition, we show 
that our data are consistent with first principles numeri- 



cal models suggesting that the Rashba-like spin-orbit ef- 
fect in biased TI junctions derives directly from the sym- 
metry inversion of topological order, and is instrumental 
in maintaining the topological connectivity of the surface 
2D band structure at large interface potentials. 

The surface state of TI Bi2Se3 crystals cleaved in ul- 
tra high vacuum (UHV) is a singly degenerate Dirac cone 
that spans a bulk semiconductor band gap of A>250meV 
(Fig. 1(a)), with a Fermi level that can be moved through 
the gap by bulk dopants such as Ca. Evaporating strong 
electron donors on P-type bismuth selenide moves the 
surface chemical potential across the band gap, achiev- 
ing a conventional diode like interface voltage of up to 
Vd~0.7V in our data from a sub-nanometer surface layer 
that is permeable to photoexcited electrons [25]. Al- 
though metal overlayers are used in this study to ac- 
commodate the shallow penetration depth of photoemis- 
sion, a thicker and less metallic semiconductor overlayer 
is likely be more practical for conventional P-N device ap- 
plications, in which case a gating potential might be re- 
quired to achieve (or exceed) the interface potentials con- 
sidered here. The cleaved Se 2_ surface of Bi2Se3 is main- 
tained at low temperature (generally T<15K), which is 
critical to provide a stable and chemically neutral sub- 
strate for these adatoms, and angle resolved photoe- 
mission spectroscopy (ARPES) performed at Advanced 
Light Source beamlines 12 and 10 is used to map the elec- 
tronic states with energy resolution better than 20 meV. 
Measurements with scanning tunneling microscopy and 
low energy electron microscopy indicate that atoms such 
as Cu and Fe do not clump on the surface, but rather 
grow in layers {if! [27J . Deposition depth is presented in 
units of monolayers (ML) to approximately represent the 
number of deposited atoms per hexagonal unit cell of the 
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FIG. 1: Topological connectivity in a TI- metal P-N in- 
terface: (a) The (blue) bulk and (green) surface band struc- 
ture of Bi2Se3 is traced above ARPES measurements of band 
structure with N-type bulk doping. A diagram above the 
ARPES image shows topological quantum numbers and the 
location of the topological surface states, (b) ARPES im- 
ages are shown for cathode (Cu-, In- and Fe-) deposited P- 
type Bi2Se3, measured after saturating the surface so that the 
binding energy of surface states no longer changed with fur- 
ther deposition. The Z2 topological quantum number iso of 
interfaced materials is indicated above, (c) The surface states 
of Cu-deposited Bi2Se3 resonate at different photon energies. 
Measurements have a <lnm penetration depth into the crys- 
tal [25[ . An integrated comparison reveals the band structure 
traced at right, with points at which bands converge labeled 
in black, including new Dirac points labeled Dl and D2. The 
Fermi surface is shown as an inset. States that evolve adia- 
batically from the original upper Dirac cone (DO) are drawn 
in a darker hue. 



surface, based on rough linear assumptions outlined in 
the online Supplemental Material [28|. 

After an estimated 1-3ML of cathode deposition on 
hole doped (P-type) Bi2Se3 the Fermi level ceases to 
move as more adatoms are deposited. The topological 
surface state is dramatically changed by the nanoscale 
P-N interface, manifesting new intersecting Dirac cones 
and a gap in the surface density of states near the original 
(DO) Dirac point after the surface state binding energy 
has been moved by 0.5-0.7eV (Fig. 1(b)), corresponding 
approximately to the chemical interface potential in nu- 
merical simulations (V^). Contrasting photoemission im- 
ages measured with different incident energies reveals the 
connectivity of energy /momentum band contours traced 
in Fig. 1(c), which must be known to understand topo- 
logical order at a Tl-metal junction. Related studies of 
N-N overlayer nanoscale interfaces are limited to smaller 
interface biases (V^^0.4V) [21M24I], and chemical poten- 
tials far from the Dirac point (DO) where topological phe- 
nomena are optimized. In this present work, our main fo- 
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FIG. 2: Topological basis for an unconventional 
Rashba effect: (a) Surface deposition of copper is shown 
on bulk P-type Bi2Se3 (0.5% Ca doped as Bii.ggsCao.oosSes), 
with approximate pre-deposition electronic structure traced 
on the first panel and energies shifted to align surface band 
features between panels. The mean free path estimated from 
inverse momentum width is labeled in angstroms for surface 
conduction bands, (inset) The Fermi surface is shown for 
0.25ML of Cu deposition with spin orientations labeled by 
arrows, (b) (black) Momentum dispersion curves of the up- 
per DO Dirac cone and Dl surface state are shown with (blue) 
fit curves that are composed of (red) r=0.06A -1 Lorentzians 
separated by 'A', (c) The Fermi surface after 2.2ML Cu depo- 
sition is shown with numerically predicted spin orientations, 
(d) (left) Large Rashba splitting in the surface states of topo- 
logically trivial BiTel is compared with (right) the strikingly 
different Rashba-like effect in topologically ordered surface 
electrons. Locations at which bands intersect or qualitatively 
appear to converge are labeled with yellow circles, (e) The 
topological connection between the valence band (VB) and 
conduction band (CB) is traced with black arrows on a DFT- 
based simulation using a coarse-grained band bending poten- 
tial from Ref. El. 



cus is the non-magnetic P-N overlayer interfaces studied 
for the first time, from which we observe the topological 
connectivity of interface states which may be seen more 
clearly than in earlier studies performed with smaller in- 
terface voltages. 

The incremental evolution leading to the appearance 
of multiple new Dirac cones [29] on cathode-deposited 
P-type Bi2Se3 is shown in Fig. [2fa). The build-up of 
positively charged adatoms first donates electrons to the 
original Dirac cone of freshly cleaved Bi2Se3, and then 
gradually binds four new extrinsic surface states that 
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intersect in two new Dirac points (Dl and D2) at the 
Brillouin zone (BZ) center. Tracing these new states 
from their points of intersection in the BZ center, we ob- 
serve that they exchange partners and appear to merge 
with a different surface band at momenta larger than 
k^O.lA -1 . This type of connectivity is also proposed in 
our previous work using Fe-deposition |l8j, but experi- 
mental evidence in the present study is stronger due to 
the Cu-deposited interface achieving sharper bands that 
are visible over a larger energy range, and due to the 
absence of potential magnetic interactions. Momentum 
distribution curves of the combined D0/D1 feature show 
that it becomes narrower approaching the Fermi level. 
Fitting the two bands with r=0.06A _1 Lorentzians sug- 
gests that they are ~0.02A _1 apart at 250meV binding 
energy and ~0.015A _1 apart at 75meV binding energy 
(Fig. [2(b)). These fits only constitute an estimation, as 
the exact broadening function required to fit these bands 
is not well known due to the parabolicity of the Dl band 
and other factors discussed below. The uncertainty in 
band separation due to the choice of broadening function 
may be quite large for cut #2 [28]. The combined feature 
has a velocity of V£>=2.5±0.1eV-A, and the velocities of 
the two bands appear to differ by just <7% over the 
0.2eV range approaching the Fermi level, in spite of the 
fact that one of the bands is induced by a surface poten- 
tial and the other is the seemingly unrelated topological 
insulator upper Dirac cone. 

The resulting Fermi surface after heavy deposition 
(Fig. [2(c)) is made up of three rings, of which the outer 
two are doubly degenerate within intrinsic width reso- 
lution. Doubly degenerate (or near-degenerate) Fermi 
surface rings cannot carry a non-zero Berry's phase [H|, 
meaning that the Berry's phase of 0#=7r by which a 
topological insulator interface is identified is obtained 
from the innermost singly degenerate Fermi surface ring. 
The outermost ring is found at k~ 1.7-2.0 A -1 , and is 
hexagonally warped due to the crystal symmetry [111 ] . 
We refer to all of these surface bands as topological sur- 
face (or interface) states, because they manifest at the 
surface of the material and all five Dirac bands must be 
observed to measure the O j b=7t Berry's phase of a topo- 
logical interface. This terminology does not preclude the 
accuracy of other terms (e.g. 2D electron gas, 2DEG), 
and the role of extrinsic surface states in topological sur- 
face band conductivity must be evaluated on a case-by- 
case basis for different interfaces as discussed below. 

As compared to the conventional Rashba states which 
can be adiabatically removed from the bulk band-gap, 
the TI Rashba states observed here are inherently re- 
quired by the nontrivial topology of the TI system. As 
shown in Fig. 2(e), in first principles simulations, a sub- 
stantial P-N junction voltage will make the original sur- 
face states no longer connect the bulk valence band to the 
bulk conduction band. Therefore, within band-bending 
models the TI Rashba-like states have to appear in order 



to make a zig-zag linkage of states that connect across the 
bulk bandgap, which satisfies the topology. Unlike the 
true Rashba effect illustrated in Fig. 2 (d, left), the spin 
splitting found in our models of Bi2Se3 is caused by the 
bulk symmetries that make the crystal a TI, which occur 
in a small symmetry inversion region at k<0.lA _1 in the 
BZ center. We note that the conventional Rashba effect 
causes weak splitting of the bands at all momenta and 
forbids exact double degeneracy away from the Kramers 
high symmetry points, but is too weak to be observed 
by eye within the plot. Therefore, interface states are 
nearly-doubly-degenerate throughout most of the BZ, as 
seen from the thick black line at k^>0.lA _1 in Fig. 2(e). 

Topological insulators are among the very few mate- 
rial systems with band structures that can be manipu- 
lated to achieve transport defined by a Dirac point, and 
most proposed applications for TIs derive from the be- 
havior of low density spin- helical Dirac surface electrons. 
In Fig. [3(c) we present three realizations of P-N over- 
layer nanoscale interfaces in which the Bi2Se3 surface 
state Dirac point is brought to the Fermi level where 
it will shape transport properties. The surface chemi- 
cal potential of an N-type sample (as-grown Bi2Se3 with 
[Se] 2+ defects) is lowered by depositing NO2 - , and a P- 
type sample (Bii.ggsCao. 005^63) is modified through de- 
position of Cu + and In + . Our data show that lowering 
the surface chemical potential of an N-type sample to the 
Dirac point has the effect of dramatically increasing the 
Dirac velocity from v jD =1.55±0.2eV-A to 2.30±0.2eV-A, 
observed from the lower Dirac cone band slope traced 
in Fig. [3(a). The velocity of lower Dirac cone electrons 
changes rapidly after deposition begins, but stabilizes af- 
ter the surface chemical potential falls below the bulk 
conduction band minimum 0.17±0.02eV above the Dirac 
point (Fig. [3(b)). The Dirac point energy is estimated 
from the point at which the bands are narrowest along 
the momentum axis, and the large error bar of ±0.2eV-A 
is assigned from the deviation between slope estimates 
over lOOmeV and 300meV windows beneath the Dirac 
point. With the Fermi level inside the band gap, bulk 
conduction electrons are strongly repelled from the sur- 
face and are not observed within lOOmeV of the Fermi 
level in the ARPES image. 

When the chemical potential lies inside the bulk elec- 
tronic gap as measured by surface-sensitive photoemis- 
sion, the lower Dirac cone velocity has the same value of 
2.3±0.2eV-A for both N-type and P-type samples, which 
is strikingly inconsistent with expectations based on the 
conventional Rashba effect. The Hamiltonian for a con- 
ventional Rashba effect induced by an external poten- 
tial is given as Hr = aji(k x z) • <r, where & repre- 
sents the Pauli matrices and olr oc Ez is the Rashba 
parameter which is proportional to the z-axis electric 
field that breaks parity symmetry. This Hamiltonian has 
a nearly identical form to the linear term defining the 
Bi2Se3 Dirac velocity, Hti = vo(k x z) • 3 [111 ], where 
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vd is the Dirac velocity. The lowest order combined ef- 
fect of conventional Rashba interactions and topological 
surface physics is therefore to change the Dirac velocity 
to a new value of V 



D 



\vd + c * Ez\, where 'c' is a con- 
stant. From this, one would expect the Dirac velocity 
change in opposite directions when the surface state is 
exposed to adatoms with opposite charge. The fact that 
the surface velocity does not show this sort of depen- 
dence is consistent with our numerical result in Fig. 2(e) 
that the conventional Rashba effect is actually very weak 
in surface deposited Bi2Se3, and suggests that other fac- 
tors such as bulk carrier density near the surface may be 
significant in creating the band instability. 

Because the lower Dirac cone velocities are roughly 
identical for all of these Dirac point interfaces, it is con- 
venient to go further in comparing their spectral line 
shapes. As shown in Fig. Efd), the energy dispersion 
curve (EDC) in the Brillouin zone center is similar for all 
of the photoemission measurements with the Dirac point 
at the Fermi level. These curves show a soft tail extend- 
ing to large binding energy, a broad peak at 40-50meV 
binding energy, and a much shallower slope approaching 
the Fermi level than would be expected from the Fermi 
Function convoluted with Gaussian resolution broaden- 
ing (blue curve). There is no simple way to obtain such an 
energy contour with only Fermi function and energy-axis 
broadening effects on the Dirac node. Because the en- 
ergies of bulk electronic states are biased significantly in 
opposite directions for Pyj-N ( e -g- Cu or In deposition) 
and Ntj-P (NO2 deposition) interfaces, it is unlikely that 
resonance with bulk electronic states is responsible for 
any of these features that are shared between all of the 
curves. 

To understand where this lineshape may come from, 
we briefly consider the asymmetrical roles of momentum- 
space and energy-axis intrinsic widths in defining the 
ARPES spectrum from 2D Dirac bands. Applying 
momentum broadening to an ideal 2D Dirac cone re- 
shapes spectral intensity in the BZ center to 1(E) oc 

\E I 

j-TT- — J 2 7^ /0 x 2 where 'E' is energy relative to the Dirac 
point and is the Lorentzian momentum broadening 
term. We have chosen momentum and energy broaden- 
ing parameters of r/e=0.0275A _1 and r#=0.0i2eV such 
that the resultant broadening vector is approximately 
normal to the Dirac cone in the 2+1 dimensional mo- 
mentum/energy space and corresponds with the intrin- 
sic width prior to broadening from experimental reso- 
lution. The characteristic EDC lineshape that emerges 
from momentum broadening is not particularly sensi- 
tive to the detailed values of and r#, so long as 
Tk is taken to be the dominant component determin- 
ing linewidths. The resulting EDC spectral function has 
a peak ~40meV away from the Dirac node (Fig. E^e)), 
and closely matches the experimental result when the 
T=15K Fermi function and weak energy- axis broadening 
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FIG. 3: Universal Dirac point kinetics and lineshape: 

(a) Deposition of anode NO 2 on N-type as- grown Bi2Se3, with 
the original slope of the lower Dirac cone overlaid on each 
panel. An inset shows the Dirac node Fermi surface after 2L 
deposition, (b) The lower Dirac cone band velocity is plot- 
ted as a function of chemical potential above the Dirac point 
as NO2 is deposited, (c) Photoemission images of bulk P- 
type and N-type samples with chemical potential set at the 
Dirac point via surface deposition of negatively and positively 
charged molecules, respectively (NO2, Cu and In). The lower 
Dirac cone velocity of pristine pre-deposition as-grown Bi2Se3 
is traced in white, (d) ARPES intensity as a function of bind- 
ing energy is plotted at the T-point for Dirac point interfaces 
shown above. A constant density of states curve cut off by the 
Fermi function is shown in blue, resembling poly crystalline 
gold, and a red curve shows simulated intensity obtained by 
broadening an ideal VD=2.3eV-A Dirac cone with Lorentzian 
parameters r fc =0.0275A _1 and r£=0.012eV. (e) Small panels 
illustrate the occupied electronic states of Bi2Se3 and poly- 
crystalline Au as a function of binding energy, demonstrating 
that momentum broadening can create a peak in intensity at 
finite binding energy. 



are also applied (red curve in Fig. [3fd)). 

The observation of a line shape representing strong 
momentum-space broadening is consistent with recent 
studies which suggest that microscopic surface domains 
on topological insulators have offset momentum symme- 
tries [31], causing them to be broadened in momentum 
space rather than along the energy axis when one per- 
forms a spatially averaging measurement such as the pho- 
toemission measurements discussed here (beam spot area 



A>200x200/im 2 ). Broadening in momentum space may 
also be understood as following from the same justifica- 
tion as broadening in energy, because electrons have fi- 
nite spatial coherence lengths, and momentum (e.g. 'k^') 
and spatial variables (e.g. 'x') are conjugate to one an- 
other. Typical bands approaching the Fermi level resem- 
ble pseudo-Voigt functions, cut by the step-like Fermi 
Dirac function, with intrinsic widths in energy and mo- 
mentum (TicFe) that are related by the group velocity v g 
as Tk=T e/v q . In most cases, there is no easy way to dis- 
tinguish between broadening that is 'applied' along the 
energy or momentum axis, as the resulting distributions 
of spectral intensity are nearly identical. 

However, close to a Dirac point singularity, applying a 
broadening convolution along the energy axis causes very 
different lineshapes than broadening bands along the mo- 
mentum axis, and the relation Tk=TE/vg no longer has 
a rigorous physical justification. States with significant 
nanometer-scale spatial 'texture' such as those observed 
by STM in Ref. [3l| provide one physical motivation for 
asymmetrical broadening, as the lifetime of these states 
can in principle be quite long (small Ye) even when 
their ensemble distribution in momentum space is very 
broad (large r^). It is well known that the momentum- 
integrated electronic density of states goes to zero at a 
Dirac point in the absence of self energy broadening. Mo- 
mentum broadening (1^) does not effect the distribution 
of electronic states along the energy axis, and thus al- 
lows the density of states at the Dirac point to remain 
zero. Energy-axis broadening (r#, imaginary self energy) 
causes the momentum-integrated density of states at the 
Dirac point energy to become non-zero, and to approach 
the Dirac point energy in a parabolic rather than linear 
fashion. 

Interfaces between conventional semiconductors and 
metals can be understood through the well known phe- 
nomenon of band bending, in which the conduction and 
valence bands of each material bend to line up at the in- 
terface. In the case of a TI interfaced with a normal 
metal, the inversion of symmetries between the topo- 
logical insulator valence and conduction bands makes it 
impossible for the bulk band structures to connect in 
such a simple way [2]. To understand the unusual inter- 
face properties observed in our data, we have created a 
Green's function implementation of the experimentally- 
based kp model in Ref. [3, to numerically simulate 
junctions at the surface of a semi- infinite topologically or- 
dered Bi2Se3 slab based on experimentally measured bulk 
electron kinetics from Ref. 17, 32| . Using a screening po- 
tential that matches the observed change in topological 
surface state electron energies in the Bi2Se3-cathode P- 
N interface (A /i =0.7eV), the simulation shows that the 
upper Dirac cone surface state electrons are still found 
in the first quintuple atomic layer of Bi2Se3, but the cor- 
responding lower Dirac cone electrons are absent, resem- 
bling a band gap as do our data (Fig. l(c, right)). 
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FIG. 4: Combining electrostatics and electronic topol- 
ogy: A numerical simulation of copper-interfaced Bi2Se3 
shows the partial density of low energy surface states at 
depths of 0-0. 5nm and 2-3nm inside the crystal. (white 
dashed lines) Electrons from the upper Dirac cone of Bi2Se3 
remain in the outermost 0.5nm, and electrons in the lower 
Dirac cone are located l-2nm deeper in the crystal. 



Looking two nanometers deeper inside the simulated 
crystal, we see that the missing in-gap electrons are part 
of a blurred continuum of relatively three-dimensional 
(3D) states derived from the bulk valence band. (Fig. 
[4j). Because these electrons are deeper in the crystal 
than the ~lnm photoemission penetration depth, pho- 
toemission measurements will show a void or "gap" near 
the DO Dirac point. This phenomenon has been modeled 
previously in Ref. [33|, in which it is determined that a 
continuum of 3D states overlapping in energy and mo- 
mentum with a TI surface state will cause the surface 
state dispersion to be reshaped with a gap-like appear- 
ance. Even though the electrons modeled in the right 
and left panels of Fig. [4] are physically displaced by sev- 
eral nanometers, there is sufficient wavefunction overlap 
between them to allow visible hybridization. No true 
3D band gap is present [33[, however the experimental 
identification of such a gap-like surface feature in our 
data is significant because a similar overlapping of states 
is expected in many theoretically predicted topologically 
ordered metals such as half-Heusler compounds [34] . Re- 
ducing the band bending potential by 50% in our model 
causes the gap-like feature to disappear [28j, consistent 
with recent measurements on limited-bias nanoscale N- 
N interfaces [22| and the absence of a gap after 0.25ML 
deposition (A /i =0.33eV) in Fig. Efa). 

Taken collectively, the P-N interface electronic struc- 
tures reported here reveal experimental ways in which the 
interaction of bulk and surface electrons reshapes surface 
band structure, and are consistent with numerical sim- 
ulations in which the symmetry inversion of topological 
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order generates a phenomenon that strongly resembles 
convential Rashba spin-splitting. Our experiment fur- 
ther suggests a characteristic lineshape of Dirac point 
interfaces, and how it can be understood in terms of the 
free path length (momentum broadening) of surface state 
electrons. 
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Supplemental Material for: Electron 
behavior in topological insulator based 
P-N overlayer interfaces 



SM I. Experimental Materials and Methods 

Spin-integrated angle-resolved photoemission spec- 
troscopy (ARPES) measurements were performed with 
35 to 65 eV photons on beamlines 10 and 12 
at the Advanced Light Source, both endstations 
being equipped with a Scienta hemispherical elec- 
tron analyzer (see VG Scienta manufacturer website 
(http:// www.vgscienta.com/) for instrument specifica- 
tions). The typical energy and momentum resolution 
were better than 20 meV and 1% of the surface BZ, 
respectively. Samples were cleaved in situ and typi- 
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cally measured at T<15K, in a vacuum maintained below 
8xl0 -11 Torr. Samples deposited with In were studied 
at a higher temperature range of T=80 to 100K. High 
quality single crystals of Ca x Bi2- x Se3 were grown by a 
process of two-step melting, as described in Ref. [l[. 

Copper, indium and iron atoms were deposited using 
an e-beam heated evaporator at a rate of approximately 
0.4A/minute for Cu and In and 0.1 A/minute for Fe. A 
quartz micro-balance supplied by Leybold-Inficon with 
sub- Angstrom sensitivity was used to calibrate the de- 
position flow rate prior to each experiment, leading to 
the possibility of some (~20%) inconsistency in calibra- 
tion between separate samples. Spectra of the adsorbed 
surfaces presented in the main paper were taken within 
minutes of opening the photon shutter to minimize po- 
tential photon induced charge transfer and desorption 
effects. No change was observed in the Cu-deposited sur- 
face energetics upon deliberately heating to 100K and 
exposing the sample to the photon beam for 30 minutes, 
suggesting that the surface was stable under experimen- 
tal conditions. 

Due to the very light coverage used, monolayer deposi- 
tion depth is defined in terms of the number of adatoms 
per surface unit cell of Bi2Se3, with a large in-plane 
hexagonal lattice constant of approximately 4. 15 A. The 
area of a surface unit cell is (4.14A) 2 x v / 3/2=14.8A 2 , 
and crystalline Fe and Cu have nearly identical densi- 
ties of 85 atoms/nm 3 , giving a naive conversion of 85 
atoms/nm 3 x0.148nm 2 = 12.6 ML/nm for each nanome- 
ter of calibrated deposition depth on the quartz balance. 
This value refers to the ratio of monolayer coverage on 
Bi2Se3 to the depth after an equal deposition time on 
the quartz nanometer in units of nanometers. Further, 
based on the change in surface energies as Cu is initially 
deposited, we estimate the sticking coefficient of copper 
on the chemically stable selenide surface to be Ss e ~~, 
similar to the low temperature sticking coefficient for Cu 
on a comparably textured oxide surface [2[ . Significantly 
larger sticking coefficients are incompatible with the ob- 
served shift in band energies near the crystal surface un- 
der light deposition, assuming that copper resides on the 
Se surface and has +1 valence. Normal-incidence cop- 
per growth on the chemically reactive surface of a quartz 
microbalance is generally assumed to have a sticking co- 
efficient of unity, giving a final calibrated deposition rate 
of l/(Ss e xl2.6 ML/nm)-2.4A/ML. This calibration is 
consistent with the decay of photoemission intensity from 
the outermost surface state bands as Cu is added. If we 
assume an ARPES penetration depth into the Cu layer of 
roughly 0.6nm at hz^=65eV incident photon energy and 
lnm at 35eV [3[, the decay of intensity 200meV above the 
original surface Dirac point is consistent with a thickness 
of 0.10±0.02 nm/ML for the first monolayer of deposi- 
tion, which is reasonable for a single atomic layer of Cu. 

For thicker layers such as the nominal "2.2ML" mea- 
surements, coverage is likely greater, because the stick- 



ing coefficient will change as a Cu layer forms. If we 
assume that the sticking coefficient changes linearly with 
coverage and reaches the sticking coefficient of the Cu- 
coated quartz microbalance after a number density of 
atoms equivalent to 2 Angstroms of bulk copper has been 
deposited, the actual coverage labeled as "2.2ML" in the 
main text will be a number of Cu atoms equivalent to 
2.35 ML of close packed [111] Cu planes (counting 0.2 
nm/ML for bulk Cu). 

Adsorption of NO2 molecules on Ca :E Bi2-^Se3 was 
achieved by controlled exposures to NO 2 gas (Matheson, 
99.5%). The adsorption effects were studied under static 
flow mode by exposing the cleaved sample surface to the 
gas for a certain time then taking data after the cham- 
ber was pumped down to the base pressure. As with Cu, 
In and Fe deposition measurements, spectra of the NO2 
adsorbed surfaces were taken within minutes of opening 
the photon shutter to minimize photon exposure related 
effects. 



SM II. Topological connectivity of surface state 
band structure 

Topological insulator surfaces can be identified by an 
odd number of surface states within the bulk band gap 
between certain high symmetry points in momentum 
space. In the case of Bi2Se3, an odd number of sur- 
face states must exist on any path through momentum 
space that connects the T- and M-points within the band 
gap. Diagrammatic examples of topological insulator and 
topologically trivial band structures are shown in Fig. ISll 
Through surface adatom deposition (Cu, In and Fe), we 
have observed the number of surface states crossing the 
Fermi level between T and M changing from one to three 
to five, preserving this odd number rule (Fig. [Sli d)). 
Exceptions to the rule can be introduced if time reversal 
symmetry is broken by magnetic perturbations. Further- 
more, surface states that are made degenerate with bulk 
electronic states (e.g. via band bending) are not explic- 
itly topologically protected [4(, as demonstrated by the 
parasitic metal band gap discussed in the main text. 

One consequence of this topological definition is that, 
if the original upper Dirac cone (DO) is shifted down- 
ward in energy so that it no longer connects to the con- 
duction band (defined by energetics deep in the crystal) 
as in Fig. [ST^d, right), the DO upper Dirac cone must 
merge with another topological surface (interface) state 
as has been observed in our data (Fig. 1(c) of the main 
text). The extrinsically bound surface states introduced 
by surface perturbations also intersect the momentum 
space Kramers points in Dirac points, and fit the topo- 
logical definition of topological surface states [4j . In Fig. 
[Sli d) and Fig. 2(e) of the main text, an electrostatically 
perturbed 12-layer GGA slab calculation (details in Ref. 
|5|) is used to show connectivity of extrinsically bound 
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FIG. SI: Topological insulator surface band connectivity: (a-b) Fermi surface diagrams are shown for a topological 
insulator and a non-topological insulator material, (c) A band structure diagram illustrates that at a topologically trivial 
surface, there are always an even number of surface state crossings between the high symmetry T- and M-points within the 
bulk band gap (here two), (d) Band structure diagrams show the GGA-predicted surface states for (left) vacuum-suspended 
Bi2Se3, and (middle, right) Bi2Se3 with positive charge on the surface. Details of the calculation can be found in Ref. 
Experimental Fermi surfaces for cathode-deposited Bi2Se3 are shown as insets, (right) The GGA simulation shows that 
topological insulator surface states become doubly degenerate (thick black line) at large momentum. 



surface states with the original Bi2Se3 Dirac cone across 
the full Brillouin zone. 



SM III. Intrinsic band broadening at a Dirac 
singularity 



In counting the number of surface states, it is impor- 
tant to know if each state is singly degenerate (one state) 
or doubly degenerate, counting as two. If the slope of 
a surface state approaching the Brillouin zone center is 
non-zero forming a Dirac point, then the lattice sym- 
metries of Bi2Se3 dictate that that state must be singly 
degenerate [6[. Empirically, we find that singly degen- 
erate states, including the original surface Dirac cone of 
pristinely vacuum-cleaved Bi 2 Se3, become doubly degen- 
erate at large momenta (e.g. Fig. ISTl d) and Fig. 1(c) in 
the main text), causing the outermost surface state Fermi 
surface rings to become doubly degenerate after lengthy 
cathode deposition (including Cu, In, Fe). For Bi2Se3 
the topological insulator state comes about as a result of 
band structure effects (termed "topological symmetry in- 
version") that occur in the Brillouin zone center, and the 
double degeneracy of surface states far from the topolog- 
ical inversion point likely represents the typical behavior 
of a topologically trivial surface-perturbed material. 



The broadening of electronic states along momentum 
and energy axes is a well known phenomenon that must 
be treated in an unusual way at Dirac singularities. En- 
ergy broadening occurs because energy and time are con- 
jugate variables, and the finite lifetime of electronic states 
due to scattering causes their width along the energy 
axis to be non-zero. By the same token, momentum 
(e.g. £ ka;') and spatial variables (e.g. 'x') are conjugate 
to one another, and states are broadened along the mo- 
mentum axes because scattering and localization restrict 
the translational symmetry of single-particle electronic 
states. For a band with small curvature, energy and mo- 
mentum broadening cannot be readily distinguished and 
it is sufficient to broaden only in energy in comparing be- 
tween theoretical band structure calculations and data. 
Dirac point singularities are a special case for which en- 
ergy and momentum broadening cause distinctly recog- 
nizable effects. As an example, it is well known that the 
momentum-integrated electronic density of states goes to 
zero at a Dirac point in the absence of self energy broad- 
ening. Momentum broadening does not effect the dis- 
tribution of electronic states along the energy axis, and 
thus allows the density of states at the Dirac point to 
remain zero. Energy-axis broadening (imaginary self en- 
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ergy) causes the momentum-integrated density of states 
at the Dirac point energy to become non-zero, and to ap- 
proach the Dirac point energy in a parabolic rather than 
linear fashion. Large momentum broadening can create 
energy dispersion curve (EDC) lineshapes that strongly 
resemble a band gap, as presented in the main text. 

As an example of how the effect of momentum broad- 
ening can be calculated, we have chosen momentum 
and energy broadening parameters of r/ c =0.0275A _1 and 
r#=0.012eV such that the resultant broadening vector 
is approximately normal to the Dirac cone in the 2+1 
dimensional momentum/energy space and corresponds 
with the intrinsic width prior to broadening from exper- 
imental resolution. The fitted lineshape is then obtained 
by broadening a Dirac cone with velocity V£>=2.3eV- 
A through Lorentzian convolution in momentum and 

energy (convoluting each point [k x o,l%o 9 Eo] with / ex 

1 



{k x -k x0 ) 2 Hky-k y0 ) 2 H^k/2) 



2 x (E-E Q y\(v E /2y)- This 'in- 
trinsically' broadened electron distribution is then trun- 
cated by the 15K Fermi function and finally Gaussian en- 
ergy broadened by the ~20meV experimental resolution 
to produce the red curve in Fig. 3(d) of the main text. 
The characteristic lineshape that emerges from momen- 
tum broadening near the Dirac point is not particularly 
sensitive to the detailed values of and r#, so long as 
Tk is taken to be the dominant component determining 
linewidths. 

The choice of broadening function is also relevant for 
other parts of the photoemission spectrum, such as the 
band intersection fitted in Fig. 2(b) of the main text. 
Fig. 2(b) uses Lorentzian fitting functions, which does 
not account for effects from the detailed momentum 
dependence of the Dl band, such as the interplay of 
energy— axis broadening with Rashba-like band parabol- 
icity. A scenario in which energy broadening dominates 
for Dl would result in a very different fitting function for 
cut #2 (see Fig. IS2]h and identifying the precise broad- 
ening function of Dl and D2 states may be an interesting 
subject for future studies. 



SM IV. Surface electron kinetics of Bi2Se3 in 
presence of external potential using k.p theory 

The k.p Hamiltonian terms from Ref. 0, 0] are 
adapted to include real space hopping parameters nor- 
mal to the crystal surface, so as to allow simulation of 
electrostatic band bending. The reduced Hilbert space 
provided by an experimentally-based perturbative model 
has the advantage of allowing a crystal to be simulated 
up to arbitrary thickness, and of capturing the correct 
bulk kinetics, which are renormalized with respect to 
DFT-based calculations by a factor of ~2 [8]. Unlike 
band bending models adopted from traditional materials 
such as GaAs [9], the modified k.p-derived model pre- 
sented here accurately reproduces the surface state Dirac 



cone and topological effects associated with the topolog- 
ical symmetry inversion. However, because the model 
is defined perturbatively around the point of symmetry 
inversion in the 3D BZ center, it does not reproduce M- 
shaped dispersions in the valence band associated with 
k z ~7r. Based on the symmetry analysis in Ref. [fiiQ, the 
effective Hamiltonian for Bi2Se3 near the Brillouin zone 
center for a quintuple layer is given by 



V 



k 2 /mi 
d + k 2 /rri2 





d + k 2 /rri2 
k 2 /mi 


ivk x + vk y 



xvkx vky 


k 2 /mi 
d + k 2 /m 2 





—ivk x + vky 
d + k 2 /m 2 
k 2 /mi 



where v is the Fermi velocity, mi and m 2 are orbital 
mass terms and the parameter d induces a band gap by 
allowing spin-independent hopping between spatial states 
in the top and bottom halves of the quintuple layer. The 
physical context of this low-order perturbative model is 
that electrons near the Fermi level in Bi2Se3 primarily 
occupy p z orbitals of Se and Bi, with parity at the T- 
point defined by the relative phase (0 or n) for electrons 
on the top and bottom halves of a quintuple-layer. The 
perturbative Hamiltonian therefore explicitly considers 
spatial states on the top (columns 1 and 3) and bottom 
(columns 2 and 4) halves of a quintuple layer, that are 
differentiated by SU(2) spin (spin up in columns 1 and 
2). The hopping between two adjacent quintuple layers 
is realized by 



/tzi \ 

t z tzl 

t zl 

V o o t z t zl J 



(S2) 



where t z is the hopping parameter and t z i breaks sym- 
metry between z-axis dynamics of the conduction and 
valence bands. To illustrate the single-Dirac-cone topo- 
logical surface states, we choose the parameters to be : 
mi = 0.125 eV^-A 2 , m 2 = -0.04 eV^-A 2 , d = -0.22 
eV, v = 2.5 eV-A, t z = 0.37 eV, t zl = -0.045 eV, match- 
ing the experimental band gap and conduction band 
mass. 

Fig. IS3l a) displays the result of our calculation of spec- 
tral weight from the first 1/2 layer (0.5nm), solved for 
this Hamiltonian. The Dirac-cone like surface states are 
clearly seen inside the bulk gap. 

Fig. IS3lf b) displays the result of a calculation of spec- 
tral weight for Bi2Se3 with a Cu film deposited on the 
surface, described by surface potentials defined in the 
next section (Section SM V). It is also well known that 
the ARPES intensity decays with the distance from the 
surface. Keeping this in mind, we combine the contribu- 
tions from first few layers to approximately match the ex- 
pected experimental penetration depth of <lnm [3]. The 
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FIG. S2: Self energy broadening and the Dl lineshape: (a) Panel 2(b,left) from the main text is reproduced, showing 
energies at which the combined D0/D1 feature has been examined along the radial momentum axis, (b) The lineshape that 
would be expected for Dl from (gray) energy-axis broadening of the parabolic band is compared with (black) a fixed width 
Lorentzian (r=0.06A -1 ). The parabolicity of Dl is approximated by fitting the black guide to the eye traced in panel (a). 
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FIG. S3: (a) Spectral weight from the top 1/2 quintuple- 
layer (top 0.5nm of the simulated crystal) without extra sur- 
face potentials. Binding energies are shifted to resemble the 
chemical potential of N-Type superconducting Cuo.i2Bi2Se3 
(b) Weighted intensity from the first few layers with extra 
surface potentials to simulate Cu deposited ARPES intensity 
for P-type Bii.995Cao.oosSe3. 



intensity displayed in Fig. IS3f b) is a combination of 50% 
contribution from the surface 1/2 quintuple-layer (top 
0.5nm), 25% contribution from the second layer, 12.5% 
contribution from the third layer and 12.5% contribution 
from the fourth layer. 



SM V. Estimating the surface potential after heavy 
cathode deposition 



Summing over the Luttinger count of all visible sur- 
face states shows that the total number of surface band 
electrons after a nominal deposition of 2 Cu atoms per 
surface unit cell on P-type Bi2Se3 is still only ~0.1e-/BZ, 
giving a good estimate of the total number of electrons 
present near the surface with symmetries derived from 
the bulk conduction band. The nominal bulk carrier 
(hole) density from Hall measurements is approximately 
10- 19 e+/cm 3 =0.002e+/BZ for the 3D Brillouin zone, 
meaning that negatively charged dopants will screen ap- 
proximately 0.002 surface electrons per quintuple layer 
(~0.95nm) of depth. The depth of the depletion region 
can therefore be approximated by the ratio of the surface 
charge density (pcu++Pss) and the bulk dopant density. 
( Pc n + nno - G / •> assuming a large bulk band mass, which 

v 0.002e /nm ' o & i 

is accurate for the valence band). 

Although Cu adatoms act as charge donors, they are 
not expected to approach a pure +1 ionization state af- 
ter the chemical potential has shifted by several hun- 
dred meV. For the case of a single fully ionized mono- 
layer of Cu (1 Cu + per surface BZ) the carrier depletion 
zone would have a depth of (l-0.1)/0.002=450nm, which 
is completely non-physical given the strong electric field 
that will be present at the surface (0.74 V/nm, even ig- 
noring changes in the Bi2Se3 dielectric constant). 

To reconcile the observed shift in surface band ener- 
gies (~0.7eV) with the screening mechanisms at play, it 
is necessary that the total surface charge of Cu atoms be 
only <0.2e + per surface unit cell. However the kind of 
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spectrum that appears in response to such a perturbation 
(middle panel in Fig. [S4l b)) is qualitatively incompatible 
with the measured ARPES spectrum in three respects. 
The points of incompatibility are: 1. The number of new 
bands derived from the valence and conduction bands 
is incorrect. 2. Spacing between bands at the surface 
is far too small, and no gap is expected to appear at 
the 'DO' Dirac point. 3. Dirac velocities of new bands 
(Rashba splitting) are far too small. All of these effects 
indicate that a typical screening model (e.g. MTFA or 
Poisson-Schrodinger methods) with no short range cor- 
rection severely underestimates the slope of the poten- 
tial near the interface and how strongly new (extrinsic) 
surface states are confined at the TI surface when the 
density of surface adatoms is large. To more closely de- 
scribe the surface kinetics, it is necessary to increase the 
slope of the surface potential in the top quintuple layers, 
to represent the short range effects such as plasticity in 
the dielectric constant and the near-neighbor interaction 
with Cu atoms. The modeled surface potential we have 
chosen, including polynomial long range screening and 
local perturbations, is given as 

V = V + Viz - V 2 z 2 + V nn (z), (S3) 

where the depth 'z' is represented by an integer index 
of the top 34 half-quintuple layers (16nm), with Vo=- 
0.35 eV, Vi =0.017 eV, and V 2 =0.0002 eV. The value of 
V2 is fixed by the bulk carrier density through the Pos- 
sion equation, Vo is determined by the change in surface 
energies, and Vi is obtained by minimizing the electro- 
static energy. The near- neighbor surface energy 'V nn ' is 
an additional perturbation of -0.41, -0.28, -0.19, -0.15 eV 
respectively for the first four 1/2-quintuple layers (2nm). 
The principle origin of this sharply sloped surface poten- 
tial in the top two nanometers derives from basic junc- 
tion physics, and is a universal source of uncertainty in 
band bending models. The dielectric constant is not well 
defined near an interface, and changes by a factor of ~20- 
40 from a strong cathode (e.g. Cu) to a semiconductor, 
meaning that the slope of the potential near the surface 
should be up to 40 times stronger than it would be with 
the dielectric constant of just the semiconductor. Our 
parameters attribute approximately 50% of the potential 
to this strongly sloped region. 

A recent study that uses a more basic quantum model 
neglecting topological order has applied a global correc- 



tion to the dipole constant (changing e r from 113 to 70) 
to partially account for this effect [9[. However a global 
adjustment of that sort is relatively coarse and misses the 
key element that the potential gradient is only strongly 
enhanced in extremely close proximity to the interface. 
Furthermore, the origin of the large Rashba-like spin 
splitting observed in our model (and present in our data 
and first principles DFT-based simulations) is the sym- 
metry inversion underlying topological order, and will be 
missed in any model that does not treat Bi2Se3 as a topo- 
logical insulator. 

This analysis of our data and recent literature suggests 
that the surface state dynamics after >1ML of deposi- 
tion (e.g. new bands and Rashba splitting) are almost 
entirely determined by near-neighbor interactions at the 
interface rather than long-range Coulomb interactions. 
Bulk-like screening with a roughly unchanged Bi2Se3 di- 
electric constant (e r ~100) is important at depths beyond 
2-3nm into the crystal, and appears to contribute approx- 
imately half of the total binding energy of the surface 
states, even though the z-axis potential gradient of bulk- 
like screening is too small to greatly effect the shape of 
surface bands beyond a rigid energy shift. When the 
surface potential is steep enough, hybridization between 
these bulk-like states at d>2nm inside the crystal and 
shallower surface bands can create a gap-like feature in 
the surface band structure, consistent with the theoreti- 
cal proposal in Ref. (Tl| . 
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FIG. S4: The importance of near neighbor surface interactions: (a) ARPES measurements of interface states on the 
surface of P-type Bi2Se3 are shown from Fig. 1(c) of the main text, (b) States within the top quintuple layer (^lnm) of Cu- 
deposited Bi2Se3 are modeled (left) using a strong near- neighbor surface interaction as described in the text, (center) assuming 
screening with only long-range Coulomb interactions and the experimental dielectric constant of e r =113, and (right) using the 
model described in the text with potentials reduced by 50% (interface voltage drop Vd^0.35V). Reducing the surface potential 
by 50% causes it to resemble the largest surface biases generally achieved across N-N interfaces, and accurately reproduces the 
recent experimental observation that N-N overlayer interfaces with biases up to Vd~0.35V manifest energy-axis broadening of 
the lower DO Dirac cone band but no gap at the Dirac point [lol ]. 



